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The single channel Anderson impurity model is a standard 
model for the description of magnetic impurities in metallic 
systems. Usually, the bandwidth represents the largest en- 
ergy scale of the problem. In this paper, we analyze the limit 
of a narrow band, which is relevant for the Mott-Hubbard 
transition in infinite dimensions. For the symmetric model 
we discuss two different effects, (i) The impurity contribution 
to the density of states at the Fermi surface always turns out 
to be negative in such systems. This leads to a new crossover 
in the thermodynamic quantities that we investigate using 
the numerical renormalization group, (ii) Using the Lanc- 
zos method, we calculate the impurity spectral function and 
demonstrate the breakdown of the skeleton expansion on an 
intermediate energy scale. Luttinger's theorem, as an exam- 
ple of the local Fermi liquid property of the model, is shown 
to still be valid. 
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I. INTRODUCTION 

The single channel Anderson impurity modelB is com- 
monly used to describe the physics of a magnetic im- 
purity in a conducting host and exhibits different types 
of behaviour such as mixed valency or the Kondo effect. 
Usually the case of a flat band is of interest, where the 
conduction electron bandwidth represents the largest en- 
ergy scale in the problem. For this particular .case the 
model has been solved using the Bethe ansatza, which 
enables one to calculate static properties, e.g. the impu- 
rity susceptibility or the specific heat. A numerical solu- 
tion of the problem has been given by Krishna-Murthy 
et all. 

In the context of the Hubbard model in high 
dimensionscl ("dynamical mean field theory"), the 
physics of the correlated lattice system is described by a 
single Anderson impurity interacting with a bathr-whose 
properties are determined in a self-consistent wayo. The 
bath is characterized by a hybridization strength of the 
order of the bandwidth and a non-constant density of 
states. For this case no exact results exist and one ei- 
ther has to rely on analytical approximations (which are, 
however, not available in the whole parameter space) or 
numerical means. 

In this paper we numerically solve the Anderson im- 
purity model in the extreme limit where the bandwidth 
is much smaller than the hybridization (narrow band sys- 
tems). In the context of the Mott-Hubbard transition in 
high dimensions such models are self-consistently gener- 



ated in one possible transition scenarioO. By analyzing 
the narrow band limit (without the additional complica- 
tions caused by the d — > oo self-consistency condition) we 
will establish two characteristic new features of Ander- 
son impurity models that to some extent will also hold 
for intermediate situations: i) The impurity density of 
states p d (u>) and the change in the density of states of 
the total system due to the impurity Ap(u>) show a very 
different behaviour since the conduction electrons react 
to the presence of the impurity. This is in contrast to the 
flat band case where always p d {io) — Ap(uj). It leads to 
an interesting crossover in the thermodynamic impurity 
properties of the system when the interaction is turned 
on. ii) The skeleton expansion, which plays an impor- 
tant role in deriving properties of the interacting system, 
breaks down on intermediate energy scales. This has im- 
portant implications for the Mott-Hubbard transition in 
d — oo as will be explained below. 



II. DENSITY OF STATES 



The Hamiltonian of the symmetric model is given by 
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with the conduction electrons Ck and the impurity orbital 
d. In the following, we will restrict ourselves to scalar 
k's, implying - if necessary - a reduction to s-waves in 
the conduction band. The impurity density of states is 
defined by 



-ImG dd {uj + ) , 
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where G dd is the retarded zero temperature Green's func- 
tion of the impurity orbital. On the other hand, the total 
change of the density of states due to the introduction of 
the impurity into the conduction band is given by 



Ap(u) 
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G kk {uj+) + G dd {u + )^ 
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where G^ refers to the Green's function of the conduc- 
tion electrons without the impurity. Using the equations 
of motion for Gkk' and G d k one finds 
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Ap(to) = — Im { G dd (u + ) 1 - ^- , 
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Therefore the total change in the DOS can be expressed 
as a function of the impurity Green's function. We in- 
troduce the notation 



E 



A(w) - iA(w) , 
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where A(w) is the hybridization function and the real 
part is given by the principal value integral A(w) = 
V J — Now we can express Ap(w) in terms of the 

impurity density of states and write 



Ap(w) = 1 - 
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We define the narrow band limit by the property 



0A(w) 
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leading to a negative coefficient multiplying the first term 
in (JsJ) at the Fermi energy. 

We will study one exemplary realization of a nar- 
row band system with a constant hybridization function 
A{u) — A and conduction band energies extending from 
—D to D (we set cf = in the sequel). This leads to 
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and hence the narrow band limit is defined by A > D. 
However, the main conclusions in the following analysis 
equally hold for other realizations of narrow band systems 
(R) as well, in particular even for systems without band 
edges at all. 

We notice that for our model the second term of 
(^) vanishes inside the conduction band. Luttinger's 
theoremEm for the symmetric Anderson model, which we 
will later verify numerically also for the narrow band case, 
ensures the "pinning" of the impurity spectral function 
at the Fermi energy at its noninteracting value 
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Therefore the first term in 
tion 



gives a negative contribu- 
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at the Fermi energy. In the following we discuss the rela- 
tion of Ap{uj) to thermodynamic properties of the model. 



III. CALCULATION OF THE IMPURITY 
SUSCEPTIBILITY 

In noninteracting fermionic models, the total density 
of states at the Fermi energy determines thermodynamic 
properties like the static susceptibility \. Introducing an 
impurity into the system induces a change in \ propor- 
tional to Ap(0). For the narrow-band Anderson systems 
analyzed here one therefore expects a negative impurity 
contribution to the susceptibility. This is obviously true 
for [7 = where one obtains the usual Pauli susceptibility 
in dimensionlcss units (ps — n = g = 1) 
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In order to determine whether this holds also for the in- 
teracting case, we have calculated x using the numeri- 
cal renormalization group method with discretisation pa- 
rameter A = 2 following Krishna-Murthy et al.B. At high 
temperatures, complete diagonalization of the logarith- 
mically discretized 5-site model yields essentially exact 
results (continuous lines), as no low energy information 
is needed in that case. The results for fixed hybridization 
A = 10 and half bandwidth D = 1 are shown in figure 
At large T, the almost free orbital leads to Ximp = g^- 
Upon lowering the temperature, we find a characteristic 
dependence on the value of the interaction. For U < 2A 
we see a crossover to a negative Ximp, which is due to 
the loss of spectral weight at low frequencies as a conse- 
quence of hybridization. After a characteristic minimum 
at finite temperature the susceptibility saturates at a neg- 
ative value as T — > 0. For U > 2A, we recover the usual 
positive Ximp which is strongly enhanced at large U/ A 
by the Kondo effect defining a new exponentially small 
energy scale Tk- 
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FIG. 1. Impurity susceptibility vs. logarithmic tempera- 
ture at hybridization A = 10 and half bandwidth D = 1 for 
different interactions U. 
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In spite of these different types of behaviour (depend- 
ing on the interaction strength), the impurity contribu- 
tion to the density of states at the Fermi energy is neg- 
ative ( |To| ) for any U . We therefore conclude that only 
in the weakly interacting case Ap(eF ) yields the thermo- 
dynamic properties of the model. In the strong coupling 
regime, the susceptibility enhancement is determined by 
the many-body resonance in pd(w) which - similar to the 
wide band case - has to be interpreted as a quasiparticle 
peak with a large effective mass. The low-temperature 
behaviour of the system for strong correlations is gov- 
erned by many-particle excitations which are not con- 
tained in Ap(uj) and, in fact, completely predominate 
over the reduction in the single-particle DOS Ap(ejr) < 0. 

IV. SPECTRAL DENSITY AND SKELETON 
EXPANSION 

Next, we discuss the spectral density of narrow band 
Anderson systems. We calculate the impurity orbital 
Green's function at zero temperature using the Lanczos 
method as implemented by KrauthO. We fix the values 
of the interaction U and the hybridization A and then 
successively reduce the bandwidth, thus taking the limit 
D — > 0. For each set of parameters we calculate the spec- 
tral function pd, using an Anderson star with 11 + 1 sites. 
For not too large U/ A the spectral density is found to be 
only weakly dependent on the number of sites. 

In the limit of small bandwidth at any finite value of 
the ratio U /A, we find a three-peak structure consist- 
ing of the atomic levels at oj = ±U/2 containing almost 
the full spectral weight and a central quasiparticle peak 
of width ~ D. Apart from numerical broadening effects 
there is no spectral weight between the peaks. This gives 
rise to resonances in the imaginary part ftt-thc self-energy, 
which can be seen in the following waySll3: As a conse- 
quence of Dyson's equation, for values of lu inside the 
gap (where spectral function and hybridization vanish) 
the self-energy £(a> + ) = K(u>) — iJ(u>) is given by 

E(w+) = w - AH - - 1 , (12) 
1 [oj) — iu + 

where we have defined T(u>) = V J de p ^t\ ■ The imagi- 
nary part of the self-energy therefore has the form 

J{oj) = tt5{T(uj)) . (13) 

For the spectral density found here T(oj) has zeroes at 
energies e ~ ±a/D and this leads to (5-functions in J(uj) 
inside the gap as showix.above. 

As argued previouslyo, these resonances are incompat- 
ible with the skeleton expansion, i.e. the self-consistent 
perturbation theory in U. Within this expansion, the full 
propagator Gdd(oj) is inserted into every diagram con- 
tributing to J(uj). As Gdd{oj) possesses spectral weight 



only on the small energy scale D <C VD, there is no possi- 
bility to generate the resonances at ±y/~D. We therefore 
conclude that in our model already at small but finite 
bandwidth D the skeleton expansion breaks down at en- 
ergies of the order y/D (strictly speaking ~ V AD, but 
note that here we have chosen A = 1 dimensionless). 

As a measure of the convergence of the expansion at 
lower energies we take the Fermi liquid properties of the 
Anderson model, especially the "pinning" of the density 
of states at its noninteracting value @. This is equiv- 
alent to the vanishing of the imaginary part of the self 
energy at the Fermi level. This property has been proven 
for a general class of systems by LuttingerEl using the 
skeleton expansion to all orders. A proof for the (flat 
band) Anderson model within unrenormalizcd-|perturba- 
tion theory was given by Yamada and Yosidaliil. In order 
to verify the pinning we focus on the central peak of 
Pd and compare different, not too strong interactions at 
constant broadening. The pinning of Pd(0) in the limit 
D — ► is evident from figure |2|. At even larger interac- 
tions (not shown here) , a narrow Kondo resonance devel- 
ops inside the band, which cannot be resolved well on a 
small cluster. From the above, though, we do not expect 
deviations from Luttinger's theorem for any interaction. 
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FIG. 2. Central peak in the impurity density of states for 
different interactions U = 0.2 and 4.0 (inset), hybridiza- 
tion A = 1.0 and bandwidth D — 10 -4 . We have used a 
Lorentzian broadening ofb=1.5xl0~ 5 at the quasiparticle 
resonance and b — 5 x 10 -8 for the poles at higher frequen- 
cies. The atomic levels at ±U /2 carrying most of the spectral 
weight are not shown. At the Fermi edge, npd{0) = 1.00±0.02 
for both values of U. 

As can also be seen in figure ^[ the "quasiparticle" reso- 
nance in Pd(oj) has an internal structure itself (including 
bound states outside the band for small U/A). In the 
limit D — > this can be described by a scaling function 

p d (u>) = f(u>/D) , (14) 

where / is independent of D. The <5-peaks in figure ||, 
dominant at small U, can be understood by approximat- 
ing the Hamiltonian as a "zero bandwidth" Anderson 
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modelEl where the impurity couples to one single orbital 
and the hybridization is given by a (5-function carrying 
the integrated weight J deA(e). The effect of the contin- 
uous hybridization band is to generate spectral density 
close to a; = and to create sidebands also visible in 
figure For larger values of U the weight of the bound 
states decreases (from the zero bandwidth model we ex- 
pect a decrease ~ 1/U 2 ) and they numerically merge 
with the continuum at small u>. 

For other hybridization functions in the narrow band 
limit (Q) the main conclusions of the above analysis re- 
main unaffected. A detailed numerical study is how- 
ever difficult due to the limited number of orbitals that 
can be taken into account using exact diagonalization. 
For example in a "narrow" band with band tails ex- 
tending to ±00 but |p| . ^> 1, the bound states for 
U = become sharp resonances. Likewise the imaginary 
part of the self-energy then contains sharp resonances at 
±0(V r D) instead of (5-functions. Still the resonances con- 
tain the same spectral weight as these (5-functions and the 
breakdown of the skeleton expansion occurs in the same 
manners. 



V. CONCLUSION 

In this paper we have studied the Anderson impurity 
model in the narrow band limit (0) using numerical meth- 
ods. We have found two new features as compared to the 
usual wide band limit, i) One observes a crossover in the 
impurity contribution to the susceptibility: For small in- 
teractions the impurity reduces the total susceptibility, 
for large interactions the impurity increases it. In fact, 
the same behaviour is also found for the specific heat, 
though this has not been discussed explicitly here. This 
crossover is in contrast to the observation that the im- 
purity contribution to the total density of states at the 
Fermi level is always negative. This quantity does there- 
fore not determine the thermodynamic properties of the 
system at large interactions, ii) Holding U and A fixed, 
the skeleton expansion breaks down for sufficiently small 
(but still finite) bandwidth D. The breakdown occurs at 
energies of order y/D and larger, while for smaller en- 
ergies no problems can be found. This shows that this 
specific breakdown of the skeleton expansion is a eeneric 
feature of narrow band Anderson impurity systemsEij. On 
the other hand, the skeleton expansion is an essential 
tool for deriving the locality of the self-energy in thc^d*- 
namical mean field theory in the Fermi liquid phaseQO. 
Its convergence also provides a sufficient condition for 
the analytic continuation to the noninteracting Hubbard 
model. At present, the question of the correct description 
of the Mott. Hubbard transition in large dimensions is un- 
der debateBBB: While NRG simulationsEJH at T = 
seem to indicate a "preformed gap" , no coexistence be- 
tween metallic and insulating solutions is found in-finite 
temperature Quantum Monte Carlo calculations^. In 



any case, the analysis presented above shows that for the 
preformed gap scenario (which naturally leads to an ef- 
fective action governed by a narrow band system in the 
sense of (0) ) one has to address the question whether such 
systems can be related to the original Hubbard model in 
large dimensions. 
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